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ABSTRACT 

X-ray reflection signatures are observed around multiple classes of accreting compact objects. Modelling 
these features yield important constraints on the physics of accretion disks, motivating the development of X- 
ray reflection models appropriate for a variety of systems and illumination conditions. Here, constant density 
ionized X-ray reflection models are presented for a disk irradiated with a very hard power-law X-ray spectrum 
(F < 1) and a variable high-energy cutoff. These models are then applied to the Suzaku data of the accreting 
X-ray pulsar LMC X-4, where very good fits are obtained with a highly ionized reflector responsible for both 
the broad Fe Ka line and the soft excess. The ionized reflector shows strong evidence for significant Doppler 
broadening and is redshifted by ^ 10'' km s^^. These features indicate that the reflecting material is asso- 
ciated with the complex dynamics occurring at the inner region of the magnetically-truncated accretion disk. 
Thus, reflection studies of X-ray pulsar spectra may give important insights into the accretion physics at the 
magnetospheric radius. 

Subject headings: accretion, accretion disks — pulsars: general — pulsars: individual (LMC X-4) — stars: 
neutron — X-rays: binaries 



1. INTRODUCTION 

Many astrophysical systems contain a hard X-ray source 
in close proximity to dense, relatively cold gas. In 
many cases the X-rays illuminate the dense gas, which 
scatters and absorbs the energy before ultimately emit- 
ting a spectrum that is imprinted with a variety of diag- 
nostic spectral features, most notab ly the Fe Kg fluores- 
cence line (e.g. , Lightman & Whi tj 119881: iGuilbert & ReesI 
[19881 iGeorge & Fabian .1991.: Matt etalJ 119911) . Such an 
X-ray reflection spectrum has been observed from accre- 
tion disks around numerous neutron stars (e.g., 'Cackett et aU 
I2OIO) and black holes of all sizes (e.g.. Miller 2007^ 
In these systems, modeling the observed X-ray reflec- 
tion features can lead to i mpo rtant constraints on the 
ionization stage jRoss et al.l 11999 : Ross & Fabian 2005j; 
iGarcfa & Kallmanl2010[). metaUicit v (Baflantvne et al. 2002), 
and densitv strucfiire feaUantvne et al. 2001, 2004, 2005) 
of the inner accretion disk. In addition, the reflecting re- 
gion of the disk may be subject to rel ativistic effects that 
will sculpt the emission features (e.g., iFabian et al.l 119891 : 
iLaoii 1199 ll) leading to measurements of disk radii and 
black hole spin (Brenne man & Revnoldsll2006l: iMilleii 120071: 
[Reynolds & Fabianii2008i) . It is clear that modeling X-ray re- 
flection features is a very important method in understanding 
the physics of accreting compact objects: however, such fun- 
damental measurements are only possible by comparing the 
observed spectra to grids of models that are computed for the 
appropriate radiative conditions of the system under consider- 
ation. 

Here, we calculate ionized X-ray reflection models appro- 
priate for accretion powered X-ray pulsars. X-ray pulsars 
are a class of X-ray binaries in which the neutron star has 
a strong magnetic field, so that accreting matter follows the 
field lines and falls into the magnetic poles, generating pulses 
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as the magnetic poles rot ate in and out of our line-of-sight 
(e.g., Nagase 1989: Bilds ten et al1ll997h . The mner accre- 
tion disk is generally believed to be truncated at the magne- 
tospheric or Alfven radius, where the energy de nsity of the 
system is dominated by the magnetic field (e.g., iFrank et al] 
120021) . Reflection of hard X-rays from the neutron star by 
the inner disk has been widely invoked to explain the ubiq- 
uitous soft excess component observed in X-ray pulsars (see 
iHic kox et al. 2004 and references therein) as well as fluores 
cent iron emission lines at w 6.4-6.7 keV (e . g.. IChoi et al 
Woo et a ll 1 1 995b fdal Fium e et al.l 1 1 9981: feurderi et al 



1994 



200C; 



2003 



La Bai-beraeLalJ|2001t |Endo et al. 2002: Naik & Paul 
2004) . Despite the ubiquity of these apparent reflec- 



tion features in observations of X-ray pulsars, there has been 
little work aimed at fitting the spectra with a self-consistent, 
physically-motivated reflection model, which is the goal of 
the present study. The next section briefly describes the cal- 
culation of the grid of reflection models that can be used for 
X-ray pulsar data, which i s then used in Se ction [3] to fit the 
Suzaku data of LMC X-4 (Hung et al 1 l2010l) . A summary of 
our conclusions is found in SectionlH 

2. CALCULATIONS 

Descriptions of the reflection calculations have been previ- 
ously pubHshed (e.g.. lRoss & Fabia n' 1993: " Ross et alj|1999l : 
,Ross & Fabian 2005), so only the main features of the model 
will be presented here. A one-dimensional, constant density 
slab of gas with hydrogen number density nn is illuminated 
by a X-ray continuum with flux Fx defined between 1 eV 
and 100 keV. Thus, for a given spectral shape, models can be 
described by an ionization parameter 



(1) 



Hydrogen and helium are assumed to be fully ionized (i.e., 
rie = 1.2nH), and the trea ted metals (C, N, O, Mg, Si, an d 
Fe) have the abundances of [Morrison & McCammonI (Il983h . 
The slab has a constant Thomson depth of tt = 15 ensuring 
that all X-ray photons interact several times even at the largest 
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values of ^. The calculation finishes, and outputs the angle- 
averaged reflection spectrum, when all 110 zones of the slab 
reach thermal and ionization balance. The ionization stages of 
the metals included in the calculation are C V-VII, N VI-VIII, 
O v-ix, Mg ix-xiii. Si xi-xv and Fe xvi-xxvii. 

The continua of X-ray pulsars can be described as a hard 
power-law with a high-energy cutoff: 



£;-rg(-Eout--B)/-Efoid 



E < Ecut 
E > Ecut, 



(2) 



where /e is the photon flux per energy interval, F is the pho- 
ton index, E'cut is the cutoff energy, and E'foid is the e-folding 
energy. The cutoff energies for most sources are typically 
< 30 keV (Hickox et al. 2004), much less than the > 100 keV 
cutoffs assumed in the pu blicly available power- law reflectio n 
models (e.g.. BaUantvne et all l200lt iRoss & Fabiaiil l2005h . 
Similarly, the photon indices of X-ray pulsars are F < 1, sig- 
nificantly harder than the reflection models calculated for fit- 
ting AGNs (Ross & Fabian 2005). These very hard continua 
also indicate that the reflection spectra will be insensitive to 
the 1 eV lower limit assumed in the calculation. 

A grid of 2376 models is calculated with parameters 
5 keV< Ecut < 30 keV (in steps of 5 keV), 5 keV< Efoid < 
30 keV (in steps of 5 keV), 0.5 < F < 1.5 (in steps of 0.1) 
and 1.5 < log^ < 4.0 (in steps of 0.5). The ionization pa- 
rameter is varied by fixing the density at rtn = 10^^ cm^'' 
(appropriate for r ~ 20 rg in a radiation pressure dominated 
accretion disk around a neutron star; = GM/c^, where 
M is the mass of a neutron star) and changing FxB Two re- 
flection spectra with different ^ are shown in Figure [T] The 
features in this spectra are qualitatively very similar to those 
used in modeling observations of AGNs. However, because 
of the high density in neutron star disks and the very hard il- 
luminating power-law, the log ^ — 3 model predicts a strong 
bremsstrahlung dominated soft excess at energies < 3 keV 
from the hot surface of the X-ray heated disk. The ther- 
mal bremsstrahlung continuum is modified by Comptoniza- 
tion in the hot atmosphere so that it approximately resembles 
the blackbody spectrum typically used to fit the observed soft 
excess. In addition, the low ionization model predicts sev- 
eral strong recombination lines from ionized C, N and O in 
the slab which, if measured, could constrain the abundances 
of the accreting gas. Interestingly, the soft excesses in X-ray 
pulsar spectra are ob served to be feature less in observations 
with CCD resolution (iHickoxet al.ll2004 ): therefore, if X-ray 
reflection is responsible for the soft excess, the ionization pa- 
rameter of the reflector must be large. 

3. APPLICATION TO LMC X-4 

To determine if the ionized reflection models provide a 
good description of the soft excesses observed from X-ray 
pulsars, the grid of models described above is fit to three 
phase-averaged Suzaku spectra of LMC X-4 (Hungetal] 
l2OT0h . LMC X-4 is a high-mass X-ray binary system con- 
sisting of a 1.29 Mf7) neutron star accreting from a 14.5 
MfT) 08 III companion (iKeUev et al.ll 19831; Ivan der Meer et al.1 
120071: iRawls et al] 120111) . The neutron star rotates with a 

^ The exact value of the density at the surface of a X-ray pulsar ac- 
cretion disk is unknown (and will likely be funct ions of the disk radius 
and accretion rate; e.g.. IShakura & Sunvae\1 [19731) . but is expected to be 
rtH S> 10'^'' cm~^, a region where calculations become computationally 
prohibitive. Recall, however, that the o-disk model gives the midplane or an 
average density, and the density at the surface is likely to be much lower 
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Figure 1. Reflection spectra from a constant density slab in'adiated by a 
power-law with a high-energy cutoff (Eq.|2j with F = 0.7, -Ecut = 10 keV 
and -Efold = 25 keV (dashed lines). The solid lines plot the reflection spectra 
for both a weakly ionized slab (log^ = 1.5) and a strongly ionized one 
(log§ = 3.0). The highly irradiated slab produces a strong bremsstrahlung 
dominated soft excess at energies < 3 keV. 

period of ~13.5 s, orb its its companion with a period of 
^1.4 days jWhitelll978l) . and exhibits a long-term periodic- 
ity of 30.5 ± 0.5 days that is believed to be caused by a pre- 
cessing accretion disk that periodicaUy obscures the ne utron 
star (e.g.. iLang et al.ll 1 98 U [Heemskerk & van Paradii sll 1 989t 
Naik & Paul 20(51. 

LMC X-4 was observed by the Suzaku X-ray observatory in 
2008 January-April, with exposures of w20 ks duration taken 
at three different phases in the superorbital cycle (^30 = 0.39, 
0.27, and 0.07 for Observations 1, 2, and 3, respectively). 
Details of the data re duction and spectral extraction are pre- 
sented in lHung et al.l (12010 ). Phase-averaged spectra were ex- 
tracted in the 0.5-10 keV energy range for the X-ray Imaging 
Spectrometer (XIS) CCD detector, and the 10-50 keV range 
for the PIN instrument of the Hard X-ray Detector ( HXD) 
and th e spectra fitted jointly for each observation. Hun g et alJ 
(I20T3) found that the broad-band Suzaku spectra are well de- 
scribed by a hard power-law with a high-energy cutoff, a soft 
X-ray excess, and a strong iron emission line, similar to pre- 
vious studies (e.g. . Woo et al. 1996; La Barbera et al. 2001]; 
iNaik & Paulll2003h . The best-fit models also included emis- 
sio n features a t « 0. 6 keV and « 1 keV that were identified 
bv lHung et aD(l20T(ih as emission from highly ionized O and 
Ne, respectively. An analysis of the energy-resolved pulse 
profiles showed that the phase offset between hard (> 2 keV) 
and soft (< 1 keV) pulsations varies with tim e, as has been 
observed for the similar X-ray pulsars Her X- 1 (iRamsav et 
^OOl fZ ane et al. 2004) and SMC X-1 (Neilsen et al. 200- 



Hickox & Vrtilek,200 5). As discussed in Hung et al.. (20101) . 



this behavior is consistent with the origin of the soft compo- 
nent as the illuminated inner region of the precessing accre- 
tion disk. 

In this paper we further explore the physical origin of the 
reflected emission in LMC X-4, by re-analyzing the Suzaku 
spect ra using our reflection model. XSPE(Z! v. 1 2.7 (lArnaudI 
11996*) is used to fit the Suzaku data. Errorbars on the best 
fit parameters are the 2a uncertainties for one parameter of 
interest (i.e., Ax^ = 2.71). 

3.1. Spectral Fits 
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Table 1 

LMC X-4 Spectral Fitting Results 



Parameter 


Observation 1 


Observation 2 


Observation 3 


Xl (385 d.o.f.) 1.13 
Afn (X 10^2 cm-2) (fixed) 0.057 
i^O.6-50 keV (erg cm-2 s"!) 8.6 X IQ-^O 


1.22 
0.057 
1.4 X 10"^ 


1.04 
0.057 
1.1 X 10"'^ 


Continuum 


r 

Ecut (keV) 
Bfold (keV) 


18.7 ±0.5 
13.3 ±0.5 


0.67 ±0.01 
18.0 ±0.4 
I5.6l0;i 


0.67 ±0.02 
17.9l°:« 
14.9 ±0.7 


Redshifted & 


: Doppler Broadened Ionized Reflection Spectrum 


log 5 

2 

o"6 kcV (keV) 


3.02 ±0.01 

034+" 0" 

19+0-07 
^•^^-0.05 


2.99 ±0.01 
034+0 004 

n -|o+o.06 


3.02 ±0.01 
036+0 012 
37+0.16 


Unblurred Reflection Spectrum 


log 5 (fixed) 


1.5 


1.5 


1.5 


Gaussian Emission Component 


E (keV) 
a (keV) 
EW (eV) 


0.94 ±0.01 
0.11 ±0.01 
129 


0.95 ±0.01 
0.10 ±0.01 
107 


0.94 ±0.01 
0.11 ±0.01 
99 


Relative Normalizations 


^XISl 
^PIN 


1.01 ±0.01 


0.95 ±0.01 
1.03 ±0.01 


1.01 ±0.01 

1 Ofi+003 


0.6-50 keV Flux Ratios 


Ionized Reflector/Power-Law 0.17 
Neutral Reflector/Power-Law 0.05 


0.09 
0.03 


0.11 
0.04 



T 0.1 



£0.01 



? 10-3 



% 1.4 

CO 

E 1.2 

0) 

E 1 

5 

0.8 

CO 

° 0.6 



■ ' ' ' 1 ' 1 




''' Till — '-^ H 1 1 ( — 1 — 1 — I— 

... 1 . 1 


SI 



5 10 
Energy (keV) 



20 



Note. — The Doppler broadening was modeled with the 'gsmooth' convolu- 
tion model in XSPEC with the Gaussian width cr{E) = erg kcv(E/6 kcV). 



The s pectral model used in the fitting is the same as the one 
used bv lHunp et al J ( 120101) except the blackbody and emission 
lines a re replaced by the ionized reflection model. Hung et aD 
(|20T^ found that both a narrow and broad Fe Ka line were 
required to fit the Suzaku data. Therefore, two reflectors are 
used in the new spectral model, with one subject to Doppler 
broadening using 'gsmooth'. The two reflection components 
and the primary power-law all have the same F, Ecut and 
i?toid- The ionization parameter of the unbroadened reflec- 
tion spectrum is fixed at log ^ = 1.5 so that it can account for 
the observed narrow, neutral Fe Ka line, but its normaliza- 
tion is allowed to vary. The results of all three fits are shown 
in Table [T] and show that in every case the ionized reflection 
model is an excellent description of the broadband spectrum 
of LMC X-4 (see Figure |2]l. In order to fit the centroid of the 
broad Fe Ka line at ^ 6.5 keV (Hung et al.l2010D . the ionized 
reflection spectrum (which produces a 6.7 keV line from He- 
like iron) must be redshifted by ^ 10^ km s^^. Interestingly, 
the fits require an additional Gaussian line with an equivalent 
width (EW) of - 100 eV at an energy of « 0.94 keV. The 
possible origin of this line is discussed below. 

As se en in Figure l2l and in agreement with the expecta- 
tions of iHickox et al.ni2004 . flie soft excess of LMC X-4 is 
almost entirely explained by a highly ionized reflection spec- 
trum, and is produced by the thermal bremsstrahlung emitted 
by the hot, ionized surface of the irradiat ed region. In addi- 
tion, the broadened Fe Ka line required bv lHung et alj (120101) 



Figure 2. The best-fitting model (top panel) and data-to-model ratio (lower 
panel) found for Observation 1 of LMC X-4. In the upper panel, the thick 
solid line plots the total model, while the individual components (the contin- 
uum, the reflectors and the Gaussian emission line) are indicated with dot- 
ted lines. In the bottom panel, the black/red/green data points are from the 
Suzaku/PVXlSl/PlN detectors, respectively. 

can be described by a redshifted and Doppler broadened ion- 
ized Fe Ka fine from the same reflector. The width of the ion- 
ize d Fe Kg line (~ .2 keV) is less than the one determined 
by iHung et al.l (120101) (~ 0.4 keV) because the fine is natu- 
rally broadened due to Comptonization (see Figure [T]l. The 
low-ionization reflector self-consistently e xplains the narrow 
O VIII and Fe Ka lines (iHun g et ani2010 l). However, the re- 
flection components are a small fraction of the total observed 
flux. The low-f reflector contributes only « 4% of the total 
flux for each observation, indicating that this emission region 
subtends a consistently small solid angle as seen from the pul- 
sar, and therefore originates from the outer disk, the accretion 
stream, or the companion star. In contrast, the highly ion- 
ized reflector contributes w 10% of the broadband flux, but 
this increases to 17% for Observation 1. As these observa- 
tions were taken at different super-orbital periods, as the inner 
accretion disk precess around the pulsar, then this changing 
fraction strongly implies that the ionized reflector arises from 
the inner regions of the magnetic ally- truncated accretion disk, 
consistent with the origin of the broadened Fe Ka line. More- 
over, the change in flux ratio indicates that our view of the 
inner reflecting region evolves with the super-orbital period, 
as expected for a warped precessing disk (Hung et al. 2010). 

The fits require an additional, broad emission component 
between 0.9 and 1 keV that does not appear in the reflection 
model (the best fit without this feature results in a Ax^ ~ 
+200). This emission is well modeled by a broad Gaussian 
component with a centroid energy of 0.94 keV, a — 0.1 keV, 
and EWrs 100 eV. These properti es are significantly differ- 
ent than the weak, narrow line that lHung et alJ (120101) needed 
at a similar energy and tentatively identified as arising from 
highly ionized Ne. In this case, a narrow emission line is 
strongly rejected by the fit, and thus favors an interpretation 
that the Gaussian component is modeling a smooth part of 
the continuum that is missing from the highly ionized reflec- 
tion spectrum. The energy and width of this feature are very 
similar to a Gaussian bum p in the XMM-Newto n RGS spectra 
of LMC X-4 analyzed by ISfeilsen et all (12009!) , who modeled 
the soft excess using a variety of different continua and persis- 
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tently observed an additional Gaussian component. Thus, this 
feature is unlikely to be an instrumental artifact, but a true un- 
modeled component of the continuum (a similar featu re was 
also f ound in the spectrum of Her X-1; Jimenez-GarateetaD 
l2002h . The energy range spanned by this emission fea- 
ture encompasses many transitions of Fe XVII-XXIII (e.g., 
iBrown et al.l 120021) . which, for this value of ^, are expected 
to be weak. However, there will likely be a steep ionization 
gradient on the surface of the inner accretion disk that is be- 
ing parameterized here by only a single ionization parameter 
In this scenario, a single ionized reflection spectrum may not 
be able to account for all the spectral complexity at < 1 keV, 
a reg ion which is ve ry sensitive to the ionization parameter 
(e.g.. iGarcia & KaUm an 2010). Whatever the origin of this 
feature, it only contributes 9% of the soft excess flux (mea- 
sured between 0.6 and 2 keV), as compared to 56% produced 
by the ionized reflector. Therefore, this unmodeled emission 
does not alter the conclusion that the highly ionized reflector 
provides a good description of the soft excess in LMC X-4. 

The redshift and Doppler broadening of the ionized reflect- 
ing gas indicates that it is subject to strong dynamical effects. 
One possibility is relativistic smearing from motions close to 
the neutron star (e.g., Fabian et al. 1989). However, tests with 
the 'rdblur' model show that good fits can only be obtained 
if the inner radius of the reflector was ^ 20-40 rg from the 
star. Identifying this radius with the magne tospheric radius 
implies a magnetic field strength of 10^ G jGhosh & Lamb! 
[1978, 1979; Frank et al. 2002), significantly smaller than the 
~ 10^"^ G field strength required to explain the stability of the 
LMC X-4 pulse period (Woo et al. 1996). Moreover, if the 
disk extended down to ^ 20-40 Vg the accretion luminosity 
from the disk would overwhelm the observed spectrum unless 
the accretion rate was < 0.07 of the Eddington rate, which is 
problematic given the observed luminosity of 3 x 10^^ erg s^ ^ . 
Thus, as indicated by its low covering factor, this gas must 
originate much farther out from the neutron star, perhaps from 
the inner edge of the disk as it interacts with the magneto- 
sphere. With a magnetic field strength of ^ 10^^ G, the mag- 
netospheric radius of LMC X-4 is ~ 2 x 10^ cm which would 
have an orbital velocity of ~ lO'' km s~^ dHung et al. 2010), 
in good agreement with both the measured redshift and the ve- 
locity width at 6.7 keV. This conclusion is supported by the re- 
sults of iNeilsen et al.l (l2009h who uncovered evidence in high 
resolution gratings spectra for multiple regions of high-speed 
ionized gas flows. It is clear that ionized reflection spectra 
not only explain the soft-excess of X-ray pulsars, but allows a 
probe of the complicated dynamics associated with the inter- 
action of the magnetosphere and accretion disk. 

4. CONCLUSIONS 

This Letter presents new reflection models developed for 
accreting X-ray pulsars. Thermal bremsstrahlung subject to 
Comptonization from an ionized X-ray heated disk is found 
to be a natural explanation for the soft excess observed in 
the spectra of LMC X-4, and, presumably, other X-ray pul- 
sars which exhibit Fe Ka lines and soft excesses. Fitting the 
Suzaku data of LMC X-4 nicely illustrates the potential of 
reflection modeling for elucidating the physics of these sys- 
tems. We found that the observed broad Fe Ka line could be 
explained as originating from a redshifted 10"* km s~^) 
highly ionized reflector subject to Doppler broadening. These 
velocities are consistent with those expected from the magne- 
tospheric radius of ~ 2 x 10^ cm dHung et al.ll2010l) . Thus, 



applying this grid of reflection models to both archival and 
future X-ray pulsar spectra will provide important constraints 
on the ionization state, metallicity and dynamics at the inner 
edge of magnetically-truncated accretion disks. 
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ful report. This work was supported in part by NSF award 
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